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This  paper  presents  a  theoretical  investigation  of  a  water- 
augmented  turbofan  engine,  one  in  which  large  quantities  of  sea  water 
are  injected  into  the  fan  duct  section.     Results  indicate  that  up  to 
three  or  four  times  dry  thrust  and  propulsive  efficiency  are  obtained 
depending  on  vessel  speed,  fan  pressure  ratio,  and  amount  of  water 
injected „     Optimum  water  injection  velocity  is  investigated.     Deviations 
from  thermal  and  dynamic  equilibrium  in  the  mixing  processes  are 
investigated  with  respect  to  their  effect  on  overall  performance. 
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NOMENCLATURE 

SYMBOL  MEANING 

a  Sound  speed    (ft/sec) 

A  Area    (ft2) 

BR  Bypass  ratio 

Cp  Specific  heat  at  constant  pressure  (ft-lb/slug-°R) 

FAR  Fuel-air  ratio 

g  Acceleration  of  gravity    (ft/sec2) 

H  Enthalpy    (ft-lb/slug) 

HVF  Lower  heating  value  of  fuel    (BTU/lbm) 

J  Joules  constant    (778  BTU/ft-lb) 

M  Mach  number 

MR  Mixture  ratio 

MW  Molecular  weight 

m  Mass  flow  rate    (slug/sec) 

P  Pressure    (lb/ft2) 

PP  Partial  pressure    (lb/ft2) ,   (atm) 

R  Gas  constant    (ft-lb/slug-°R) 

S  Entropy    (ft-lb/slug-°R) 

SFC  Specific  fuel  consumption    (lbfu/lb-hr) 

ST  Specific  thrust    (lbf-sec) 


lbm 


T  Temperature  (°R) 

TH  Thrust  (lb) 


SYMBOL 


V 


MEANING 


Velocity    (ft/sec) 


V  Mean  velocity    (ft/sec) 

WGR         Water  to  gas  generator  air  ratio 


X 

Specific  humidity 

2T 

Ratio  of  specific  heats 

1 

Efficiency 

P 

Density    (slug/ft3) 

SUBSCRIPT 

MEANING 

A 


Air 


B 

Burner 

C 

Compressor 

D 

Gas  generator  diffusor 

f 

Saturated  liquid 

fg 

Change  by  evaporation 

fu 

Fuel 

F 

Fan 

FD 

Fan  diffusor 

g 

Saturated  vapor 

i  Refers  to  property  at  station  (i) ,  i  =  1,2.  .  . 

I  Refers  to  state  reached  by  isentropic  process 

N  Gas  generator  nozzle 

P  Propulsive 

REFA  Reference  of  air 


SUBSCRIPT  MEANING 

S  Static 

T  Total,  Turbine 

W  Water 

V  Vapor 


CHAPTER  I 
INTRODUCTION 

With  the  present  day  emphasis  on  nuclear  propulsion  and  its  far 
reaching  effects  on  a  modern  Navy,  advances  in  other  types  of  marine 
propulsion  may  appear  less  necessary  and  certainly  less  glamorous. 
However,  the  significance  of  small  to  medium  displacement  craft  in  a 
limited  conflict  has     been  forcibly  brought  home  by  the  role  such  craft 
are  playing  in  Vietnam  and  in  other  troubled  and  far-reaching  corners 
of  the  globe.     One  of  the  major  features  of  the  limited  conflict  is  the 
speed  and  mobility  of  participating  units,  generally  with  respect  to 
land  and  air  forces;  however,  naval  forces  must  soon  follow  with  faster 
vessels.    With  existing  types  of  propulsion,   speed  increases  must  be 
paid  for  with  dramatic  weight  increases  often  outweighing  the  speed 
increase.     Consequently,  one  of  the  most  recent  fields  of  endeavor  in 
marine  propulsion  has  been  to  develop  a  system  that  supplies  enough 
power  to  attain  high  speeds  and  yet  is  lightweight  by  comparison.    This 
paper  deals  with  one  of  the  most  recent  proposals,  the  water-augmented 
turbofan  engine. 

Gas  turbines  are  currently  being  used  to  drive  conventional  screw 
propeller  vessels  and  more  exotic  systems  such  as  the  water-jet  driven 
hydrofoil  (PCH-2)  of  the  U.  S„   Navy,     However,  little  serious  thought 
has  been  given  to  using  turbines  for  jet  propulsive  power  since  the  high 
exit  velocities  inherent  in  gas  turbines  (even  in  turbofan  engines) 
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produce  characteristically  low  propulsive  efficiency  and  are  thus  too 
impractical.     Conventional  thrust  augmentation  by  water  injection  was 
considered;  however,  this  generally  involved  injection  of  water  into  the 
compressor  air  flow  or  into  hot  turbine  exhaust  gases.    In  the  case  of 
the  former,  high  purity  water  would  be  necessary  because  of  the  contact 
with  vital  moving  parts.    In  both  cases  the  high  temperatures  of  the  gas 
generator  would  vaporize  most  of  the  water,  and  the  energy  required  to 
accomplish  this  would  negate  the  thrust  increases  (1).    A  mixed  flow 
turbofan  where  fan  duct  and  gas  generator  exhaust  through  the  same 
nozzle  would  have  the  same  deficiency  since  high  temperatures  in  the 
combined  nozzle  would  again  vaporize  most  of  the  water.     Davison  and 
Sadowski  (2)  proposed  a  novel  method  of  injecting  sea  water  into  the 
relatively  cool  flow  of  a  fan  duct  where  the  fan  and  gas  generator  exhaust 
through  separate  nozzles.    As  a  consequence,  very  little  water  is 
vaporized,  and  no  moving  parts  are  contacted,  making  the  use  of  sea- 
water  feasible.     See  Figure  1. 

The  purpose  of  this  paper  is  to  optimize  the  design  point  parameters 
of  bypass  ratio,  water  to  gas  generator  ratio,  and  those  velocity, 
temperature,  and  pressure  ratios  that  directiy  affect  the  performance  of 
this  system  with  respect  to  specific  thrust  and  propulsive  efficiency.    In 
particular,  it  is  desired  to  determine  to  what  extent  departures  from 
thermodynamic  and  dynamic  equilibrium  in  the  mixing  processes  affect 
performance  and  what  critical  limits  exist,  if  any.    It  would  then  have  to 
be  established  experimentally  that  a  mixing  duct  with  two-phase  flow  in 
the  actual  conditions  would  satisfy  these  requirements. 
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This  is  both  necessary  and  important  for  two  reasons.     First,  the 
work  of  Davison  and  Sadowski  (2)  ,  which  constitutes  the  only  major 
effort  on  the  subject,  deals  with  one  particular  engine  and  one  fixed 
set  of  operational  parameters .    In  other  words ,  they  have  shown  it  to 
be  feasible  for  one  particular  engine  but  have  made  no  effort  to  find 
under  what  conditions  such  a  system  would  be  optimum,  and  such 
information  is  necessary  for  the  design  of  a  prototype.    Second,  an 
experimental  mixing  duct/fan  nozzle  section  should  be  constructed  to 
experimentally  verify  the  theoretical  findings  of  this  paper,  especially 
with  regard  to  acceptable  limits  for  velocity  and  temperature  lags 
between  water  and  air  and  also  mixing  duct  lengths  necessary  to  achieve 
them. 

The  data  obtained  in  the  paper  were  obtained  from  a  computer 
program  designed  to  analyze  the  performance  of  the  water-augmented 
turbofan  engine  while  parametrically  varying  the  critical  design  parameters 
so  that  graphical  and  tabular  displays  of  thrust  and  efficiency  variation 
could  be  obtained  and  optimum  points  determined. 

Briefly,  the  remainder  of  the  thesis  will  present  a  description  of 
the  model  used,  develop  the  equations  upon  which  the  computer  program 
is  based,  present  methods  of  analyzing  the  data,  and  conclude  with 
results,  discussion,  conclusions,  and  suggestions  for  further  work  to 
be  done  in  this  area. 
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CHAPTER  II 
DESCRIPTION  OF  MODEL  USED 

2.1  Physical  Layout 

Figure  2  is  a  block-schematic  diagram  of  the  operating  system 
depicted  in  Figure  I.    Important  points  to  consider  are  that  the  fan  duct 
system  is  completely  separate  from  the  gas  generator  system  with 
separate  exhausts  and  that  water  injection  occurs  in  the  mixing  duct 
entrance  after  the  fan,  thus  contacting  no  moving  parts.    All  stations 
are  numbered  in  Figure  2;  and  when  a  numerical  subscript  is  involved 
in  an  equation,  the  quantity  thus  subscripted  is  measured  at  that  station 
in  the  engine.    It  should  be  noted  that  the  figures  do  not  show  fuel  tanks, 
water  injection  pumping  systems,  or  any  other  components  contained  in 
the  actual  engine. 

2.2  Assumptions 

Because  of  the  large  number  of  variables  in  an  analysis  of  this 
type,  many  simplifying  assumptions  had  to  be  made  to  limit  computational 
time.    The  air  was  assumed  to  be  behaving  as  a  perfect  gas  with 
constant  specific  heats.    As  is  common  practice  in  rough  design  work 
in  turbomachinery  applications,  Cp  was  assumed  constant  at  one  value 
for  low  temperatures  and  constant  at  a  second  value  for  high  temperatures. 
Gamma  was  then  computed  from  the  appropriate  value  of  Cp.     For  the 
purposes  of  this  investigation,  a  plot  of  Cp  versus  temperature  of  the 
gas  was  made  from  data  in  reference  (3).    The  results  of  this  were  two 
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values  for  Cp  for  the  two  temperature  conditions.     For  the  combustion 
chamber,  which  is  the  transition  between  the  two  extremes,  an  arithmetic 
mean  of  the  two  values  of  Cp  was  used.    This  was  justified  on  the  basis 
that  the  plot  of  Cp  versus  temperature  was  roughly  a  straight  line  in  the 
region  of  interest.    Viscous  effects  were  ignored  with  the  exception  of 
the  influence  of  the  water  on  the  air  in  the  fan  duct.    It  was  assumed 
that  no  heat  was  transferred  other  than  in  the  combustion  chamber,  and 
no  work  was  transferred  other  than  in  the  turbine,  compressor,  and  fan. 
Mixing  duct  area  was  assumed  constant,  and  it  was  also  assumed  that 
the  engine  would  have  only  converging  nozzles.    In  other  words,   both 
gas  generator  and  fan  exit  at  atmospheric  pressure;  however,  if  this 
produced  an  exit  Mach  number  greater  than  one,  the  values  were 
recomputed  to  exit  at  a  higher  pressure  with  sonic  velocity.    The  last 
major  assumption  made  in  the  initial  analysis  was  that  the  air  entering 
the  system  was  dry  air  with  no  relative  humidity  effects.    For  ease  in 
calculations  one  slug  of  air  was  assumed  to  be  entering  the  gas 
generator  per  unit  time.    As  a  consequence,  bypass  ratio  was  also  the 
mass  flow  through  the  fan  section  per  unit  time,  and  an  increase  in 
bypass  ratio  keeping  other  parameters  constant  was  an  increase  in  fan 
section  size  holding  gas  generator  size  constant. 
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CHAPTER  III 

ANALYSIS  OF  DRY  TURBOFAN 

With  zero  water  injection,  the  analysis  of  the  system  is  developed 
from  H-S  diagram  relationships  and  basic  equations  of  engineering 
science.     This  chapter  and  the  equations  contained  in  it  are  arranged  in 
sections  similar  to  those  contained  in  the  computer  program  for  the  dry 
turbofan  engine. 

3.1   Gas  Generator  Diff user 

Ambient  total  pressure  is  determined  from 


tf-i  ..M*-1 


Vro-    -Psol*    TMo]  (1) 

so  that  defining  diffuser  efficiency  as 

1=     ^        -  <D__  <2> 


it  follows  that 


TV.  •     ?so    +    \\ft*  -  ^o) 


Since  no  work  is  done  in  the  diffuser, 


and  enthalpy  is 


Ur,  =  CPc  TTI-  I*,    +    H« 


r-z, 


where  HREpA  is  the  enthalpy  of  the  air  at  TsQ. 
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(3) 


PT|     =      Tro    =         Tso    I    l  +  ^^o]  (4) 


(5) 


3.2   Compressor 


The  ideal  (isentropic)  temperature  at  compressor  exit  for  a  specified 
pressure  ratio    — is 


\> 


ri 


Tm  =    TT* 


\*TX 


Pn 


(6) 


Define  compressor  efficiency  as  the  ratio  of  isentropic  work  input  to 
actual  work  input;  thus 


H 


Ti 


Cpcl    T-T2.1      ~      \T\)    I      Ut 


(7) 


T2. 


1  VAtz.     -       M  V36FA 


Cpc 


\  so 


3.3   Combustor  (Burner) 

Recall  from  Section  2.2 


2..0 


Pe> 


(8) 


TT^  is  fixed  by  maximum  allowable  turbine  inlet  temperature 
Therefore 


Ht 


Cf0      Ttj    -     »T1        +      Hti 


Writing  the  energy  equation  across  the  burner, 


(9) 


If  it  is  assumed  that  injection  temperature  of  the  fuel  is  equal  to  air 
temperature  at  that  point,  and  the  energy  needed  to  vaporize  the  fuel 
is  included  by  choosing  the  lower  heating  value  of  the  fuel,  then 
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FAK   ~- 


C?H  [    \r$   -  TtZ  ) 

V  HVF.fr- J     -    C?w(T-r5-TT2.] 


(10) 


3.4  Turbine 

For  the  entire  system  work  input  must  equal  work  output;  therefore , 
HT4  may  be  calculated  from 


and 


Hti"   Wr 


■*•    en* 


Hti-  V4T(e 


Tr4--      W3 


-r  _       (Ht3-JAta    ) 


(11) 


(12) 


£.Pw 


Turbine  efficiency  is  the  ratio  of  actual  work  output  to  isentropic  work 


output,   so  that 


\TA-I    "       W3 


Ht5     -     vAt^) 


(13) 


An  expansion  whether  isentropic  or  actual  is  to  the  same  pressure; 


therefore , 


IPt*     -      ^tai      - 


?T3 


I  T4I 
TT3 


*-» 


(14) 


3.5    Gas  Generator  Nozzle 

Since  no  work  is  done  in  the  nozzle, 

Tt  A-    ~       ^ T  ^> 

Knowing  that  the  flow  exits  to  ambient  pressure  fixes  the  pressure 
?S5 


ratio 


Thus 


V- 


T4- 


SSI 


T4 


'5  5 


(15) 
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Define  a  nozzle  efficiency 


^  TV*    -  TSE>x  (16) 


Equation  (16)  may  be  solved  forTqq  and  exhaust  velocity  determined 
from 


v*5  =yic?c(TTr~  ls\) 


(17) 


Sound  speed  is  _  /.     o  ~r   ~ ' 

Therefore  exit  Mach  number 

KBi;     ^  as) 

Note  the  comments  in  Section  2.2  concerning  procedure  when  IVL  !!>1 


3  .  6  Fan  Diffuser 

Specifying  a  separate  fan  diffuser  efficiency     *^J  FD   ,  the  equations 
for  the  fan  diffuser  correspond  to  equations  (2)  through  (5)  with  subscript 
6  substituted  for  subscript  1 . 

3.7  Fan 

With  a  corresponding  fan  pressure  ratio  and  fan  efficiency, 

h_    ,  the  fan  equations  are  identical  to  the  compressor  equations  of 
Section  3.2  with  the  subscripts  7  and  6  substituted  for  the  subscripts 
2  and  1  respectively. 

3.8  Fan  Mixing  Duct 

No  work  is  done  in  the  mixing;  therefore, 


/ 
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\T7      '-         We>  (19) 


?T6    = 


1  ?T7  (20) 


PT7 


3.  9  Fan  Nozzle 

The  equations  of  Section  3.5  are  valid  for  the  fan  nozzle  also  with 
subscripts  9  and  8  substituted  for  5  and  4  respectively  and  the  proper 
fan  nozzle  efficiency  specified.    It  should  be  noted  at  this  point  that  the 
entropy  at  all  stations  is  calculated  from  relationships  similar  to 

St;   =    cPt  ufej     -    ttUfej   +  S,eFA       (2i) 

where  Srefa  is  the  entropy  at  TgQ  anc*  PgO  . 

3.10    Specific  Thrust,  Propulsive  Efficiency,  and  Specific  Fuel 
Con  sumption 

Thrust  is  determined  from  the  momentum  flux  equation 
TH  -.     rnTV^    -  ^cVo   *  rnpN^-Vo)  +  J^sfes  -^°)   (22) 
Specific  thrust  based  on  the  total  mass  flow  of  air  entering  the  system 


is 


ST  =       <%  [ri 


rr\c  +  m  F 


(2  3) 


Propulsive  efficiency  is  a  measure  of  how  well  the  system  is  converting 
the  change  in  kinetic  energy  of  the  working  medium  into  thrust  power  of 
the  vehicle  and  can  be  shown  to  be 
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V 


TV*  -  V0 1 


p  TH-vo    -v    mT(\^-Vo)      -k    m^(V/v^  -  v/q)  (24) 

Specific  fuel  consumption  is  the  number  of  pounds  of  fuel  required  to 
produce  one  pound  of  thrust  continuously  for  one  hour  and  can  be  shown 
to  be 


SFC    =        FAR-    &c  i    SkOQ  ■  <y 


(25) 


TH 
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CHAPTER  IV 
FIRST  WATER  INJECTION  ANALYSIS 

The  injection  of  water  in  the  mixing  duct  does  not  affect  the  results 
of  Sections  3.1  through  3.7.     The  equations  contained  in  these  sections 
are  still  valid  for  the  case  of  water  injection.    The  new  analysis  is 
based  on  solving  the  problem  in  the  mixing  duct  and  fan  nozzle  and 
evaluating  the  new  specific  thrust  and  efficiency  from  the  information. 

In  this  first  water  injection  analysis  it  was  assumed  that  none  of 
the  water  in  the  mixing  duct  and  fan  nozzle  sections  was  vaporized. 
The  rationale  behind  this  assumption  is  that,  according  to  reference  (2), 
much  less  than  one  per  cent  of  the  water  is  vaporized. 

4 . 1    Definitions 

The  following  terms  will  be  used  throughout  the  chapters  dealing 

with  water  injection: 

Bypass  Ratio  (BR)  =      m  F  /  rv%  c 

Water  to  Gas  Generator  Ratio  (WGR)      =      rr\  w  /  rmc 
Mixture  Ratio  (MR)  =      rn  w/  /  nn  P 

The  mean  velocity  at  a  station  is  defined  after  Elliot  (4)  as  the 

velocity  that  the  air  and  water  would  attain  if  they  were  allowed  to 

reach  equilibrium  at  that  station  with  no  loss  in  momentum. 

t^^V/\      ■*"      ^vj\/w       =      [**\*   ■*■  A/^|\l  (26) 


2.4 


thus , 


Mean  Velocity  (V)    =  \/a    +     H^-Vw 


(27) 


4.2   Mixing  Duct 

The  momentum  equation  for  the  assumptions  of  the  paper  is 


0nF  r  rr\w 


cW   --    -Ad 


? 


(28) 


The  area  occupied  by  the  air  is  m  p  /^     .  . 

The  area  occupied  by  the  water  is        cow  /n      w 


Thus,  the  total  flow  area  is    A  = 


rn 


V\R 


(29) 


i 
Combining  equations  (29)  and  (28)  and  utilizing  the  fact  that  the  mixing 

duct  is  assumed  to  be  of  constant  area,  equation  (28)  can  easily  be 

integrated  and  becomes 

_  i  -     if         \         W      1 MR     \ 

V6  =     V7     +     [?s7    "     ^8|      ,  4  nRh  fWA    ^     pwVw  (30) 


The  energy  equation  is 


r^p 


^S7 


z 


V 
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Equations  (30)  and  (31)  have  five  unknowns:  PS8,  VA8 ,  VW8 ,  TS8,  and 
TW8 .    With  three  equations  lacking  for  a  solution,  it  was  necessary  to 
specificy  ratios  of  unknowns  to  serve  as  equations.    These  three 
remaining  equations  were 
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TSS8WA    =  TW&/TS6 


(32) 


VR3WA 


VWb/VAS 


(33] 


PSR87 
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(34) 


These  admirably  suited  the  eventual  task  of  determining  the  effect 
of  deviations  from  equilibrium  in  this  mixing  process  on  overall 
performance . 

4.3   Fan  Nozzle 

The  momentum  equation  cannot  now  be  integrated  to  yield  an  equation 
of  the  form  of  equation  (30)  since  variation  of  area  and  pressure  in  the 
nozzle  is  unknown.     Several  attempts  were  made  to  provide  a  second 
equation  either  from  known  facts  or  from  assumptions  concerning  area 
variation  in  the  nozzle  with  varying  results.    The  problem  was  rendered 
academic  when  the  underlying  assumption  of  no  water  vaporization  was 
found  to  yield  erroneous  results.     (See  Section  4.4.) 


4.4   Conclusion  of  First  Attempt 

A  continuity  equation  was  introduced  into  the  mixing  duct 
calculations  in  order  to  reduce  the  number  of  variable  ratios  specified, 


It  was  introduced  in  the  form  of 


h-,    =    A8 


or  more  specifically 
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The  problem  was  now  overspecified  with  more  equations  than  unknowns. 
When  two  of  the  three  ratios  in  equations  (32) -(34)  were  used  with 
continuity,  momentum,  and  energy  as  defined  in  equations  (30),   (31), 
and  (35),  ridiculous  values  were  produced  for  the  unknown  variables. 
Careful  analysis  of  the  equations  involved  revealed  that  they  were  correct 
for  the  given  assumptions  indicating  that  the  assumptions  needed 
re-evaluation.    The  net  result  was  that  the  assumption  concerning  no 
vaporization  of  the  water  was  not  correct,  and,  indeed,  vaporization  of 
the  water  was  found  to  play  a  large  role  in  determining  the  solution. 
Admittedly  a  very  small  percentage  of  water  is  vaporized,  but  for  two 
reasons  its  effect  is  significant.     First,  the  entalpy  of  water  vapor  is 
roughly  one-thousand  BTU's  per  lbm.  greater  than  that  of  liquid  water. 
Second,  the  amount  of  water  present  with  respect  to  air  (MR)  is  a  number 
ranging  in  value  from  roughly  ten  to  four  hundred.    Thus  the  contribution 
in  an  energy  equation  of  vaporized  water  had  to  be  taken  into  account 
along  with  the  resultant  change  in  water  mass  flow. 
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CHAPTER  V 

SECOND  WATER  INJECTION  ANALYSIS 

The  second  analysis  is  similar  to  that  of  Davison  and  Sadowski  (2) . 
However,  it  differs  in  that  their  work  involves  the  assumptions  of  thermal 
and  dynamic  equilibrium  in  the  mixing  duct,  dynamic  equilibrium  in  the 
fan  nozzle,  and  zero  pressure  change  in  the  mixing  duct  while  this 
paper  takes  into  account  possible  variations  in  these  quantities. 

5  .  1   Steam  Table  Properties 

It  was  necessary  in  the  second  analysis  to  calculate  the  properties 
of  enthalpy  of  the  liquid  (Hf) ,  enthalpy  change  by  vaporization  (Hfg) , 
enthalpy  of  the  saturated  vapor  (Hg) ,  entropy  of  the  liquid  (Sf) ,  entropy 
change  by  vaporization  (Sfg) ,  entropy  of  the  saturated  vapor  (Sg) ,  and 
partial  pressure  of  the  vapor  (PPv)  „    These  are  all  complicated  functions 
of  temperature.    A  least  squares  method  was  used  to  fit  data  points 
read  in  from  the  steam  tables  (5)  in  the  range  of  interest,  and  it  was 
found  that  a  cubic  approximation  gave  accuracy  down  to  the  last 
significant  figure  in  the  steam  tables.    The  exception  was  partial  pressure 
of  the  vapor,  the  actual  equation  of  which  was  simple  enough  to  be 
included  in  the  program  directly.    Since 

and 


S<>  -     5*    +     S^  <37> 
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it  was  necessary  only  to  compute  the  properties  Hf,  Hfg,  Sf,  and  Sfg, 
thus  determining  Hg  and  Sg.    The  equations  generated  by  the  four  least 
squares  approximations  and  the  PPv  calculation  are 
HF  =  -5.2244716xl02  +  (1 .  1719422)  (Tw)  -  (3  .  1882153xl0~4)  (Tw)2 

+  (1.9421466xl0~7)(Tw)3  (38) 


HFG  =  1.440548xl03  -  (1 .  1017804)  (Tw)  +  (1.0967588x10    3)(Tw)2 

-(7.4163358xl0~7)(Tw)3  (39) 


SF  =  -1.6947235  +  (5.325x10  3)(Tw)  -  (4.7499752x10  6)(Tw)2 

+  (1.8842757xl0~9)(Tw)3  (40) 


SFG  =8.812883  -  (2.442480x10    2)  (Tw)  +  (2.8390143x10   5)(Tw)2 

-(1.2441932xl0"8)(Tw)3  (41) 


Equations  (38)  through  (41)  are  valid  in  the  range  of  temperatures 
from  510  °R  to  660  °R.    Tw,  the  temperature  of  the  liquid,  is  in  degrees 
Rankine,  and  the  units  of  output  for  the  curve  fit  are  BTU/lbm  for  enthalpy 
and  BTU/lbm-°R  for  entropy.     In  the  equations  that  follow  it  will  be 
assumed  that  Hfg,  for  example,  has  already  been  converted  into  ft-lb/ 
slug. 
log10(PPv)  =  log10(Pc)  -    &    [A+_1X^CX3_1  (42) 


where 


Pc  =  218.167  atm 


X  =  (647.27  -  T)     °K 


T  =  (Tw  -  492)  (9/5)  +  273.16      °K 
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A  =  3.2437814 

B  =  5.86826xl0~3 

C  =  1.1702379xl0~8 

D  =  2.1878462x10" 3 
and  the  answer,  PPv,  will  be  in  atmospheres.    Again,  in  the  equations 
that  follow,  it  will  be  assumed  that  PPv  has  been  converted  into  lb/ft 

5  .  2    Mixing  Duct 

The  specific  humidity,  X,  of  the  mixture  is  defined  as  the  ratio  of 
the  mass  of  water  vapor  to  the  mass  of  dry  air  in  the  system.     From 
perfect  gas  equations  for  dry  air  and  water  vapor,  it  can  be  shown  that 

/ vWj 1/iaWv    j/Tss    \  ,,,. 

XS-(?si-PM6](    W,   ](Tve    I  (43) 

It  was  assumed  that  vapor  temperature  was  equal  to  liquid  temperature. 
Taking  into  account  vaporization,  the  energy  equation  is  now 

mF    I   \AS1    V    V^iM     +  r*^  [U^     +    Vvo^j     =     ^FL.V^)    (44) 


•v      m^sj 


The  momentum  equation  (29)  and  its  integrated  result  (30)  are  still 
valid  for  the  case  of  vaporization.    The  problem  is  now  one  of  three 
equations:  (44),    (43),  and  (30)  with  six  unknowns:  VW8 ,  VA8 ,  TW8 ,  TS8, 
PS8,  andX8.     By  again  specifying  ratios  of  unknowns,  equations  (32) 
through  (34),  a  solution  is  possible.     However,  an  explicit  solution  is 
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not  possible  since  equations  (4  3)  and  (44)  contain  involved  functions 
of  TW8  through  the  enthalpy  and  partial  pressure  terms.    VA8  and  VW8 
may  be  solved  using  equations  (30),   (33),  and  (34).    The  remaining 
equations  may  be  manipulated  in  such  a  manner  as  to  yield  a  function 
of  TS8  only  (called  FUN)  which  is  equal  to  zero.    A  Newton  Rhapson 
iterative    solution  may  then  be  applied 

TS8.=TS8,    .     -    FUN/DFUN  (45) 

where  TS8.  is  the  jt^1  approximation,  TS8.       is  the  (j-l)*1*1  approximation, 
and  DFUN  is  the  derivative  of  the  function  with  respect  to  TS8,  both 
evaluated  at  TS8.    ,  .    When  the  correct  value  of  TS8  is  reached,  FUN  will 
equal  zero,  and  TW8  and  X8  may  be  calculated  from  equations  (32)  and 
(43). 

5  .  3   Fan  Nozzle 

Along  the  same  reasoning  as  in  Section  5.2,  specific  humidity  at 
the  fan  nozzle  exit  is 

/     MwVWTsq 


Aq.       - 


wWA      Ts/s 


(4  6) 


Note  again  vapor  temperature  was  assumed  equal  to  liquid  temperature. 
With  varying  area,  the  momentum  equation  cannot  be  solved  unless 
a  geometric  configuration  is  specified.    Therefore,  the  assumption  was 
made  that  the  net  change  in  entropy  through  the  nozzle  is  zero.    That  is 

(s<\  -   s*  L      *  (s«  -  s&]       -  o      <47> 
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Since  exit  pressure  is  atmospheric,  the  equations  (48)  and  (46)  have  TS9, 
TW9,  and  X9  unknown.     By  specifying 


TSR9WA  =  TW9/TS9 


(4  9) 


and  solving  equations  (46)  and  (48)  for  a  function  of  TS9  equal  to  zero, 
the  Newton-Rhapson  iterative  method  again  yields  values  for  TS9,  TW9, 
and  X9 .    The  energy  equation  may  be  written 


™  f=      W^  8 
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Let  VR9WA  =  VW9/VA9 

and  solve  equations  (50)  and  (51)  for  VW9  and  VA9 
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5.4   Specific  Thrust  and  Propulsive  Efficiency 

Thrust,   specific  thrust,  and  propulsive  efficiency  are  defined  as 
before,   but  account  must  be  taken  of  the  varying  mass  flow  rates  in  the 
mixing  sections. 
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It  should  be  noted  that  ambient  air  velocity  was  assumed  equal  to  ambient 
water  velocity.     In  other  words,  wind  and  current  were  neglected. 


33 


CHAPTER  VI 
RESULTS  AND  DISCUSSION 

6 .  1    Organization  of  Output 

The  results  obtained  focus  on  two  basic  indicators  of  water-augmented 
turbofan  performance.    These  are  the  ratio  of  specific  thrust  for  the  water- 
augmented  engine  to  specific  thrust  for  the  dry  engine  and  the  ratio  of 
propulsive  efficiency  for  the  wet  engine  to  propulsive  efficiency  for  the 
dry  engine.     Hereafter  these  two  parameters  will  be  referred  to  as  thrust 
ratio  and  efficiency  ratio.     The  final  data  concern  two  broad  areas  of 
interest,  the  effect  of  variations  in  general  system  parameters  on  thrust 
ratio  and  efficiency  ratio  holding  mixing  section  equilibrium  parameters 
fixed  at  1.0,  and  the  effect  of  variations  in  mixing  section  equilibrium 
parameters  on  thrust  ratio  and  efficiency  ratio  holding  all  other  system 
parameters  constant. 

Subsequent  figures  show  plots  of  these  variations  and  also  the 
effect  of  water  injection  velocity  on  system  performance. 

6.2   General  System  Parameters  vs.  Thrust  and  Efficiency  Ratios 

The  parameters  originally  thought  to  have  the  most  significant  effects 
on  thrust  ratio  (TR)  and  efficiency  ratio  (  *[  R)  are  water  to  gas  generator 
air  ratio  (WGR) ,   ship  speed  (VO) ,  and  fan  pressure  ratio  (FPR),     Performance 
is  also  improved  by  increasing  the  bypass  ratio  (BR),  but  this  parameter 
was  not  thoroughly  investigated  since  its  limit  is  a  function  of  design 
considerations  only. 
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Figure  3  presents  thrust  ratio  versus  WGR  for  varying  VO  and  FPR. 
The  more  water  injected,  the  greater  is  its  influence  on  the  air  as 
evidenced  by  a  decrease  in  fan  nozzle  exit  velocity  (VA9) .     From  the 
thrust  equation  ,     equation  (6  3)  ,  it  can  be  seen  that  a  point  is  reached 
with  increasing  WGR  where  the  negative  contribution  of  the  momentum 
drag  of  the  water  outweighs  the  positive  contributions  to  thrust.    That  is, 
the  momentum  flux  into  the  system  eventually  increases  more  rapidly 
than  the  momentum  flux  out.    Thus,  in  general,  the  curves  increase, 
reach  a  maximum  or  peak  value,  and  then  fall  off  as  more  water  is  added. 
Note  that  for  a  higher  value  of  VO,  values  of  TR  are  much  lower.    This 
is  because  of  the  increase  in  inlet  momentum  drag  of  the  air  and  water 
with  respect  to  other  thrust  contributions.    Also  note  that  at  the  higher 
values  of  VO,  the  curves  peak  at  a  lower  WGR.    Again,  this  is  because 
the  inlet  momentum  drag  of  the  water  predominates  over  the  increase  in 
thrust  much  more  rapidly. 

In  comparing  Figure  3  with  the  curves  obtained  by  Davison  and 
Sadowski  (2) ,  close  correlation  is  found  in  shape  and  numerical  values 
with  differences  attributable  to  the  fan  pressure  ratio.     In  this  paper  fan 
pressure  ratio  is  defined  as  the  ratio  of  total  pressure  aft  of  the  fan  to 
total  pressure  before  the  fan,  thus  making  mixing  duct  static  pressure  a 
variable  with  ship  speed,  while  Davison  apparently  defines  fan  pressure 
ratio  as  the  ratio  of  static  pressure  aft  of  the  fan  to  ambient  pressure, 
since  his  mixing  duct  static  pressure  is  constant  for  all  ship  speeds. 
Figure  3  shows  the  fan  pressure  ratio  to  be  an  extremely  important  design 
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parameter  when  water  injection  is  considered.     The  values  of  thrust 
ratio  obtained  clearly  indicate  that  excellent  values  of  peak  thrust 
augmentation  could  be  obtained  in  the  order  of  upwards  of  three  times 
dry  thrust  to  over  one  and  one-half  times  dry  thrust. 

Figure  4  presents  efficiency  ratio  versus  the  same  independent 
parameters  with  a  constant  fan  pressure  ratio.    The  shape  of  the  curves 
obtained  is  nearly  identical  to  those  in  Figure  3  with  roughly  similar 
numerical  values.     Efficiency  ratios  peak  at  slightly  higher  values  of 
WGR  than  the  corresponding  thrust  ratio  indicating  that  efficiency  ratio 
is  slightly  less  sensitive  to  the  increase  in  momentum  drag  of  the  water. 
No  particular  peak  value  predominated  over  the  range  of  velocities  for 
either  thrust  ratio  or  efficiency  ratio.     However,  for  a  specified  velocity 
and  fan  pressure  ratio,  a  nearly  identical  value  of  optimum  WGR  was 
found  for  both  thrust  and  efficiency  ratios.     Significant  was  the  fact  that 
at  higher  velocities  this  value  was  relatively  low  indicating  that  high 
speed  craft  are  equally  if  not  more  suitable  for  water  injection  that  those 
operating  at  lower  speeds. 

6.3  Equilibrium  Parameters  vs.  Thrust  and  Efficiency  Ratios 

For  the  analysis  of  variations  in  mixing  duct  and  fan  nozzle 
equilibrium  parameters  on  thrust  and  efficiency  ratios,  craft  speed  was 
fixed  at  fifty  knois,  fan  pressure  ratio  at  1.5,  bypass  ratio  at  4,  WGR 
at  50,  and  compressor  pressure  ratio  at  13.8.    The  four  equilibrium 
parameters  chosen  were  the  ratio  of  water  to  air  temperature  at  mixing  duct 
exit  (TW8/TS8)  and  fan  nozzle  exit  (TW9/TS9)  ,  the  ratio  of  water  to  air 
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velocity  at  mixing  duct  exit  (VW8/VA8) ,  and  fan  nozzle  exit  (VW9/VA9) . 
Four  runs  were  made  holding  all  but  one  of  the  equilibrium  parameters 
fixed  at  1 .  0  and  varying  the  remaining  one  through  its  maximum  possible 
range,    (Figures  5  through  £) .     For  a  slightly  more  meaningful  output  the 
two  temperature  ratios  were  varied  simultaneously  with  the  velocity 
ratios  held  constant  (Figure  9)  and  likewise  the  velocity  ratios  varied 
together  with  temperature  ratios  fixed  (Figure  10). 

Figures  5  and  7  indicate  that  for  all  other  equilibrium  parameters 
fixed  at  unity,  equilibrium  at  the  mixing  duct  exit  is  either  detrimental 
or  has  no  effect  on  thrust  and  efficiency  ratios.    This  result  was  at  first 
thought  erroneous  but  then  rationalized  as  being  correct  (although 
unrealistic).     Note  that  for  these  curves,  equilibrium  conditions  are 
required  at  the  exit  of  the  fan  nozzle.    Thus  by  requiring  less  equilibrium 
in  the  mixing  duct,  the  mixing  is  being  forced  to  take  place  in  the 
nozzle,  where  it  occurs  isentropically,  resulting  in  higher  efficiency 
and  thrust  ratios . 

The  results  obtained  in  Figures  6  and  8  are  much  more  significant 
as  they  represent  not  only  physically  possible  states  but,  in  the  case 
of  Figure  6,  physically  probable  ones.     Figure  6  shows  that  thermal 
equilibrium  at  the  nozzle  exit  is  definitely  not  desirable.     Peak  thrust 
occurs  for  values  of  TW9/TS9  between  1.0  and  1.05,  which  can  amount 
to  as  much  as  twenty-five  or  thirty  degrees  temperature  difference.    At 
greater  values  of  this  temperature  ratio,  thrust  falls  off  again  although 
much  more  gradually,  indicating  that  a  wise  operating  condition  would 
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be  to  the  right  of  the  maximum  point.     Efficiency  ratio  is  not  affected 
to  nearly  the  extent  of  thrust  ratio  although  it  also  increases  with  a 
divergence  from  equilibrium  and  in  fact  never  does  begin  to  peak. 

Figure  9  plots  thrust  and  efficiency  ratios  versusAemperature  ratio 
at  nozzle  exit  varying  concurrently  with  temperature  ratio  at  mixing  duct 
exit.    In  this  case,  the  predominance  of  the  nozzle  exit  temperature 
ratio  is  again  felt  as  perfect  thermal  equilibrium  is  seen  to  be  the  worst 
possible  condition  with  a  peak  va lue  occurring  shortly  after  the 
equilibrium  condition.     Figure  10  shows  the  same  type  of  plot  with  both 
velocity  ratios  varying  and  both  temperature  ratios  held  constant.    It 
reflects  the  fact  that  the  velocity  ratio  in  the  mixing  duct  makes  little 
or  no  difference  on  thrust  holding  the  other  ratios  fixed  in  that  it  is 
almost  identical  to  Figure  8.    Again  departure  from  equilibrium  means 
a  decreasing  thrust  ratio  and  efficiency  ratio. 

6.4    Other  Results  of  Analysis 

One  of  the  more  significant  results  of  the  paper  was  discovered 
because  of  initial  failures  to  assume  realistic  data.    Initially  a  constant 
water  injection  velocity  was  specified  regardless  of  craft  speed.    It 
was  felt  that  this  could  be  achieved  with  proper  spray  nozzle  design  and 
a  pumping  system  to  provide  the  necessary  pressure  differential  when 
the  speed  of  the  vessel  through  the  water  did  not  create  enough  head  to 
inject  the  water  at  the  desired  velocity.    This  constant  injection  velocity 
specification  gave  absurdly  high  values  of  thrust  ratio  not  believed  to 
be  possible.    With  the  new  assumption  that  the  injection  velocity  is  a 
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fraction  of  the  vessel  speed  (obviating  the  use  of  pumps)  data  were 
obtained  which  closely  correlate  with  previous  work  (2) .     However,  the 
fact  remained  that,  for  a  higher  injection  velocity,  values  of  specific 
thrust  were  extremely  high.    Several  TR  versus  WGR  curves  for  varying 
injection  velocity  are  presented  in  Figure  11  with  equilibrium  parameters, 
bypass  ratio,  and  fan  pressure  ratio  fixed.    It  was  then  postulated  that 
if  a  pumping  system  were  incorporated  to  provide  the  higher  injection 
velocities,   such  a  system  could  be  used  in  emergencies  and/or  for 
important  missions  when  a  burst  of  additional  thrust  was  required. 
Naturally,  the  power  required  to  run  such  a  pumping  system  would  reduce 
the  propulsive  efficiency,  but  this  would  be  well  offset  by  the  otherwise 
unobtainable  increase  in  thrust.     Note  that  this  would  not  require  more 
water  injection  but  only  higher  injection  velocities.     In  effect,  it  would 
be  a  type  of  sea-borne  "afterburner. "    Such  a  device  in  intelligence 
ships  such  as  the  USS  PUEBLO  would  be  of  unquestioned  value  for  use 
in  emergencies.     Perhaps  retractable  hydrofoils  could  be  incorporated  to 
be  used  only  when  on  the  high  thrust  mode  of  operation. 
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FIGURE  11 
THRUST    RATIO  VS.  WATER  TO  GAS  GENERATOR  RATIO 
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CHAPTER  VII 
CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  STUDY 

7. 1   Conclusions 

The  proposed  concept  of  water  injection  appears  to  be  a  feasible 
means  of  increasing  thrust  and  propulsive  efficiency.     Peak  values  of 
thrust  ratio  and  efficiency  ratio  occur  when  momentum  drag  of  the 
injected  water  predominates  over  the  thrust  due  to  additional  mass  flow. 
These  peak  values  occur  at  lower  values  of  water  to  gas  generator  ratio 
for  increasing  vessel  speed  because  the  momentum  drag  of  the  water  is 
proportional  to  vessel  speed. 

Peak  values  for  thrust  ratios  occur  slightly  before  peak  values  for 
efficiency  ratios  since  the  efficiency  is  less  sensitive  to  inlet  momentum 
drag  of  the  water. 

Static  pressure  in  the  mixing  duct  plays  an  important  role  in  the 
slope  of  the  curves  and  thus  determines  the  effectiveness  of  the  water 
addition  scheme  at  various  ship  speeds. 

The  equilibrium  parameter  in  the  mixing  process  which  appears  to 

have  significance  in  affecting  system  performance  is  the  ratio  of  water 

temperature  to  air  temperature  at  the  fan  nozzle  exit.     Peak  performance 

with  regard  to  thrust  ratio  is  obtained  when  TW9/TS9  is  in  the  range 

of  1 .  02  to  1 .  05  .     Considerably  better  values  of  thrust  ratio  occur  on 

the  side  away  from  equilibrium  as  opposed  to  closer  to  equilibrium. 

Efficiency  ratio  continues  to  improve  as  TW9/TS9  moves  further  away 

from  equilibrium. 
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Normally,  injection  velocity  would  be  a  fraction  of  the  vessel 
speed  determined  by  how  well  the  ram  scoop  inlet  is  designed.    If  a 
pumping  system  is  installed  such  that  the  water  might  be  injected  at  a 
higher  velocity,   short  periods  of  high  thrust  could  be  experienced  with 
a  corresponding  decrease  in  propulsive  efficiency,   such  as  in  an 
aircraft  afterburner. 

7.2   Suggestions  for  Further  Study 

There  is  a  great  deal  of  further  work  required  before  an  experimental 
setup  can  be  constructed-,  much  of  this  work  can  be  accomplished  with 
this  program  and  modifications  to  it. 

The  ambient  conditions  were  fixed  throughout  this  analysis  and  an 
investigation  into  the  variations  of  system  performance  for  different  sea 
water  temperatures  and  air  temperatures  and  pressures  could  be  made. 
Since  vaporization  of  the  water  in  the  mixing  sections  was  of  such 
significance,  the  effect  of  relative  humidity  in  the  air  should  be 
investigated.    It  should  also  be  possible  to  choose  a  friction  factor  for 
the  mixing  duct  and  through  the  use  of  the  Fanno  flow  equations 
determine  the  pressure  drop  in  the  mixing  duct  rather  than  specify  it  as 
was  done  in  this  analysis.    The  question  of  water  injection  velocity  and 
a  pumping  system  "afterburner"  should  be  investigated  with  regard  to  the 
decreases  in  efficiency  caused  by  the  additional  power  required  to  run 
the  pump. 

Last  and  most  important,  it  remains  to  design  and  construct  an 
experimental  fan  section  that  can  investigate  the  effects  of  (a)  the 
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equilibrium  ratios,    (b)  static  pressure  in  the  mixing  duct,  and  (c)  injection 
velocities  of  the  air  and  water  on  thrust  and  efficiency  and  to  correlate 
the  experimental  results  with  the  theoretical  values  obtained  in  this 
program  and  in  this  analysis. 
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APPENDIX  I 

COMPUTER  PROGRAM  NOMENCLATURE 

Al.l       Input  Parameters 

Card  1: 

ETAD  Gas  generator  diffuser  efficiency 

ETAC  Compressor  efficiency 

ETAB  Burner  efficiency 

ETAT  Turbine  efficiency 

ETAN  Gas  generator  nozzle  efficiency 

ETAFD  Fan  duct  diffuser  efficiency 

ETAF  Fan  efficiency 


Card  2: 


DSO  Ambient  density    (slug/ft3) 

PSO  Ambient  pressure    (lb/ft2) 

TSO  Ambient  temperature    (degrees  Rankine) 

HVF  Heating  value  of  fuel    (BTU/lbm) 

TREFA  Reference  air  temperature    (4  92  deg.  R) 


Card  3: 


AMN7       Air  injection  Mach  number 
TWO  Sea  water  temperature    (deg.  R) 

PTR32        Pressure  ratio  across  burner 


Card  4: 


NVO 

Number 

NCPR 

Number 

NFPR 

Number 

NTT 

Number 

NBR 

Number 

NWGR 

Number 

NMDPR 

Number 

NMDVR 

Number 

NMDTR 

Number 

NFNVR 

Number 

NFNTR 

Number 

of  vessel  speeds 

of  compressor  pressure  ratios 

of  fan  pressure  ratios 

of  turbine  inlet  temperatures 

of  bypass  ratios 

of  water  to  gas     generator  ratios 

of  mixing  duct  pressure  ratios 

of  mixing  duct  velocity  ratios 

of  mixing  duct  temperature  ratios 

of  fan  nozzle  velocity  ratios 

of  fan  nozzle  temperature  ratios 
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Card  5: 

GC  Cool  temperature  value  of  ratio  of  specific  heats 

GH  Hot  temperature  value  of  ratio  of  specific  heats 

Card  6: 

VO(I)  Values  of  vessel  speed  (ft/sec),  I  =  l,NVO 

Card  7: 

PTR21(J)  Values  of  compressor  pressure  ratio  (lb/ft2)  ,1=1  ,NCPR 

Card  8: 

PTR76(M)  Values  of  fan  pressure  ratio  (lb/ft2)  ,  M  =  1,NFPR 

Card  9: 

TT3(K)  Values  of  turbine  inlet  temperature  (deg  R) ,  K  =  1,NTT 

Card  10: 

BR(L)  Values  of  bypass  ratio,  L  =  1,NBR 

Card  11: 

WGR(N)  Values  of  water  to  gas  generator  ratio,  N  =  1  ,NWGR 

Card  12: 

PSR87(IJ)  Values  of  mixing  duct  pressure  ratio,  IJ  =  1,NMDPR 

Card  13: 

VR8WA(IK)         Values  of  mixing  duct  velocity  ratio,  IK  =  1,NMDVR 
Card  14: 

TSR8WA(IL)      Values  of  mixing  duct  temperature  ratio,  IL  =  1,NMDTR 
Card  15: 

VR9WA(IM)       Values  of  fan  nozzle  velocity  ratio,  IM  =  1,NFNVR 
Card  16: 

TSR9WA(IN)      Values  of  fan  nozzle  temperature  ratio,  IN  =  1,NFNTR 
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AI .  2     Output  parameters 

VA8  Air  velocity  at  mixing  duct  exit  (ft/sec) 

VW8  Water  velocity  at  mixing  duct  exit  (ft/sec) 

VA9  Air  velocity  at  fan  nozzle  exit  (ft/sec) 

VW9  Water  velocity  at  fan  nozzle  exit  (ft/sec) 

TS8  Air  static  temperature  at  mixing  duct  exit  (°R) 

TW8  Water  static  temperature  at  mixing  duct  exit  (°R) 

TS9  Air  static  temperature  at  fan  nozzle  exit  (°R) 

TW9  Water  static  temperature  at  fan  nozzle  exit  (°R) 

STHRUS  Specific  thrust  (lb  thrust/lbm/sec) 

AMN5  Gas  generator  exit  Mach  number 

AMN9  Fan  nozzle  exit  Mach  number 

ETAPRO  Propulsive  efficiency 

SFC  Specific  fuel  consumption  (lb  fuel/lb  thrust-hr) 

FAR  Fuel-air  ratio  (lb  fuel/lb  air) 

X9  Specific  humidity  at  fan  nozzle  exit  (lb  vapor/lb  dry  air) 
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MAIN 
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APPENDIX  II 


COMPUTER  PROGRAM  AND  SAMPLE  OUTPUT 


THOMAS  C.  KNLDS°N 
aj^T  TURBOFAN 
SAMPLE  DATA  RUM 
?-■    MARCH,  1  Q  t.  q 

RFAL  mr 

DIMENSION  V0(6)  ,  PTR?1  (  1  j  )  ,  PTR76I  IT  )  ,TT3(  5)  ,  HP  (  1  ■'  )  ,WGR(  lr  )  ,  P<;PR7(  1 
*if  VR8WAUJ)  ,TSR8WA(  m.VRStfAtn)  ,TSR9WAt  10) 


STFAM  TABLES  FUNCTIONS  ANn  DERIVATIVES 


*4? 

HF 

-35 

SF 

*y 

SF 

OH 

OH 
OS 
TS 

RE 
FO 
RF 

FP 
RF 
FD 
RE 

RE 
Ffl 

pr 

kF 
RE 

3  F 

Kf 

RF 
RE 
RE 

RC 
RE 
RF 
FH 

R  = 

G  = 
CJ 

CP 

CP 

TW 
HH 


UX)  =  - 

1466*XX 

G  (  x  x  »  =  l 

3*XX**3 
( XX)=-1 
3188927 


19 


522.<+4716  +  1.  171<3422*XX  -  ."'.  '""*188?153*XX*XX  ♦  .  T" ' 

**3 

44,":. ^4R     -     1.1      178r'4*xy     ♦     .  "~  ir367ci3P*xx*XX     -     .PCTT"  7^163 


.694723 

5  7*--  XX** 

G(  XX  )=e.  81?8  8'4 
4 1 '-)  3  ?  *  X  X  *  *  3 

171Q42 

- 1  •  l  r  1 7 

53?5n 


F{  XX 
FG(  X 
H  XX 
F  G  (  X  X  )  = 
AD(5 
RMAT 
Al)(  5 
RMAT 
AD(  5 
RMAT 
ALU  5 
RMAT 
AD(5 
RMAT 
AQ(  5 
AIM  5 


AO( 
Ai)( 
AO( 
AO( 
AO( 
AO( 
AO( 
AD( 
AO( 

RMAT 
1717 
3^.2 
=77R.26 
H  =  r,M*R  / 

C=GC*R/ 
?  =  Twn 

FFA=  CP 


5     +     .  *  jr^2Vrr  1*XX     -     ,■   ?  ">  v-474QQ7c57)hxX*XX     +     .^A^O 
3 
-     .  '  ?44?48fH*  XX     +     .  ^  >.' ^?R3C)''  1  4^*XX*  YX     -     ,^'V""] 


=  1 
>  = 


?     - 


»637643"6*XX    +    .  T'.V*^3(  5  ^?^439«*XX*xy 


1) 
7F 
2) 
Fl 
3) 
3F 
4) 
11 
5) 
2F 
6) 
6) 
6) 
61 
6) 
6) 
M 
ft) 
61 
6) 
6) 
R 
M 


804    +     .v  P  219351  76*XX    -    .  "I  '  f  ^r  2  ??40^r  "74*XX*X  X 

)'J)     -     •  '•<<"  OC  r,Q4^Q9^"4*XX     «■     .  ^C  :'1C>'ir  C^6'S78271  *XX*XX 

4R"1  ♦  .rr.^r  <567?T?«6*XX  -  .  t  :  V  ">"*  3  7^  ?5706*X  x*  v  y 

ACtETA8tETATtFTANtETAFD«ETAF 


ETAD.ET 
5.2) 

DSOfPsn 

0.6.4F1 

4MN7, J  a 
8.2) 

NVO,NCPR»fl'FP0»vJTT.NPR,NWO,'.NM|;)PRtNMnvP,NMnTR.MFNVR,NFNTO 

13) 

GCtGH 

9.6) 

( V0<  I  )  . 

( PT^  21( 

{ PTP76( 

( TT3( K ) 

(  H  R  ( I  )  , 

( WGP( N) 

(  PSF«7( 

{  VRF'wAI 

(  T  S  P  8  ta  A 

{ VP9WA( 

( TSR9WA 


, TSnfHVF , TRFF A 

').  1  ) 

C, PTR32 


1=1, NVO) 
J),J=1,NCPR) 

v )  ,  m  =  1 f MFPR ) 

,  K  =  1 ,  N  T  T ) 

L=l , NBR) 

, N=l , NWGP ) 

IJ  )  t IJ  =  1,NMDPR ) 

IK  )  ,  I K  =  1  ,  NMDVR  ) 

(  IL  )  , I L  =  1  .NMOTR ) 

IM) , [M=l ,NFNVP ) 

(  IN)  ♦  I  «M=l  ,  NFNTR  ) 


(  GH-1 .1  ) 
(GC-l.-i  ) 

c*< tsh  - 


T  P  c  (-  \  ) 
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G  LEVEL  1,  MOD  o 


MAIN 


DATE  =  6BPR1 


14/^7/57 


SO 

7 


1 


HR 
SR 
WP 
FO 
WP 
FO 
WP 
FO 
WP 
FO 
WP 
FO 
*'E 
WR 
*ET 
11  FO 
*6, 
*T« 
*tT 
*-D 

*,2 

*^ 

WR 

Fn 

DO 
CV 

vw 

AO 
AM 

VW 


1? 


EFW  = 
EFA  = 
ITE( 
R"AT 
ITE( 
RMAT 
ITE( 
RMAT 
ITE( 
RMAT 
I  TF  ( 
RMAT 
FF  IC 
ITE( 
AFD, 
RMAT 
•FT- 
,T35 
3,  »T 
EG  R 
,/tT 
SPEC 
ITE( 
RMAT 

36 
=  VO( 
0  =  OV 
=  SOR 
Njn  =  o 
7  =  0. 


6,50 
CI' 
6,7) 
(/// 
6.  8  » 
<//T 
6,9) 
(//T 
6,10 
(/// 
IFNC 
6,1  1 
ETAF 

<//t 

LB/S 
t  'SR 
SO  = 
•  ,T  7 
74,  ■ 
IFIC 
6,12 
(//T 
/,T3 
1  =  1, 
I  )*1 


WO)*CJ*G 

) 
) 

///T4-),  •  TURBOFAN  ENGINE  CYCLF  PROGRAM') 

44, 'WITH  WATER  INJECTION*) 
UDE  =   0  FT* ) 

ENT  PRODERTIES' ,T35, 'REFERENCE  PROPER T I FS ', T74, 
A,FTAD,PSO,SREFA,FTAC ,TSO,HRFFW,FTAB,ETAT,ETAN, 


47,  •  Al.  TIT 
) 

/T3, • AMRI 
IES'  ) 

)DSO,HRFF 


3, 'DSO    = 
LUG' ,T74, 
FFA=     «,F1 
•  ,F9.4,T 
4, 'ETAR    = 
ETAFD    =     • 
ATIONS'  ) 
)AVN7, TW7 
3, 'A"N7     = 
, 'Tw7    =     • 
NVO 
.6PB944 


•  ,F9. 6, T21,  'SLUG/TO.  FT •  ,  T35 , • HRFF A  =  ',F11.1,T5 
•FTAP  =  SF6.2,/ ,T3, «PSO  =  •  , F9. 3  ,  T21 ,  •  L R/SO.  F 
1. 1 ,TS6, 'FT-LB/SLUG-DEG  R',T74,'ETAC  =  »,F6.2,/ 
21,'DEG  R  »,T35,'HREFW  =  •  ,  F  1 1 . 1 , T 56, • FT-L B/SLUG 
•  ,F6.2,/,T74, 'FTAT  =  •  ,  Ft . 7 , / , T74 , • FT  AN  =  «,F6 
,F5.2,/ ,T74, 'ETAF    =     • , F6. 2 , / , T3 . ' WAT FR     TNJECTIO 


'  ,F7.2,/,T3, 'VWO    = 
,F7.2,T19, 'DEG    R' ) 


VO 


KNOTS 


T(GC*P*TSO i 

V/AO 

R   *    vwr 


GAS    GENERATOR     DIFFUSES 

PTO       =  PSP*(  (l.O+l  (GC-l.Gl/2.0  )*AMN0*AMNO)f-*(GC/(GC-i,r  )  )  ) 

PT1=ETA0    *(PTO-PSO)     ♦     PSO 
TTl=TSO*(i.C ♦((GC-l.O»/2.0)*AMNO*4MNO) 

STl  =  CPC*AL0G(TTl/TSn)-R*AL0G(PTl/PS0)4-SRFFA 
HT1=CPC*<TT1-TS0)+MRFFA 

FAN    DIFFUSE^ 

PT6=ETAFD*( PTO-PSO)     +    PSn 

TT6=TS0*( 1 .C*( IGC-l.f ) /2. ^ |*AMNn*AMMn) 

ST6=CPC*AL0G(TT6'TSP)-R*A10G(PT6/PS0)+SRFFA 

HT6-CPC*<  TT6-TS0UHREFA 

GAS    GENERATCF     COMPRESSOR 

DO  35  J=1,NCPR 

TT2 I=TT1*PTP21 { J »**( IGC-l.fi/GCI 

0HT2ll=CPC*(  TT2I-TT1) 
HT2=DHT21 T /ETAC  ♦  HT1 
TT2=(HT2-HREFA) /CPC  +  TSO 
PT2=PTR21( J)*PTl 
ST2=CPC*ALOG(TT2/TSO)-R*ALOG(PT2/PSO) 


+  SRFFA 


GAS  GFNFRATOP  BURNER 


00    34    K=1,NTT 

PT3=PTR32*PT2 

CPR=(CPH    ♦    CPO/2.C 

HT3=CPR*( TT31K J-TT2)     ♦    HT2 

ST3=CPB*ALOG(TT3(K ) /TT2)     -    R * AL 0G( PT3 /PT2 )     ♦    SPEFA 

FAR=1.0/( ( FTAB*HVF*G*CJ)/(CPH*(TT3(K)-TT2) )-l.DI 

FAN 

DO    33    M=l ,NFPR 

TT7I=TT6*<  PTR 76 {")**<  (GC- 1.0) /GO  ) 

DHT76I=CPC*( TT7I-TT6) 

HT7=DHT76I /FTAF  +  HT6 

TT7=(HT7-HREFA  ) /CPC  *■     TSO 
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MAIN 


DATE  =  6R081 


14/37/57 


WV7 

WRI 

14  FDR 

*TIO 

*NLE 

*Y  ( 


*****#***<-**********  **  ft*******************************  ) 

=O.P*VC< I ) 

TE(6.l4)VC(I),PTR21(J> , PTR76I M ) , TT 3 ( K ) .WV7 

MAT(////Tln.,  »VO  (KNOTS)  IS  »,  F5.  1  ♦/ »T  1')  ,•  COMPRESSOR  PRESMJRF  RA 
IS  '  ,F5.  1, /,  TIC, 'FAN  PRESSURE  RATIO  IS  •  ,  F5.  1 , / , T 10 , ' TOR R I NF  I 
T  TEMPERATURE  (DEC.  Rl  IS  ' ♦ F 7. 1 ♦ / , T 1 0 , • WATFR  INJECTION  VFIHTTT 
KNOTS)  IS  «  ,F5.  1, /////) 


IS 


GAS  GENERATOR  TURRIMF 


DO  11 
RY=RR 
CMF=1 
F"F  =  B 
TVF  =  C 
DHT43 
HT^=H 
DHT43 
HT4I  = 
TT^I  = 
TTA  =  T 
PT4I  = 
PT4  =  P 
IF(  (  P 
ST4  =  C 


L  =  l 
(L) 

•  - 

Y 

Mp     ♦ 
=  (CM 

T3-D 

I  =  OH 

HT3- 

TT3( 

T3(K 

PJ3* 

T4I 

T4-P 

PH*A 


NBP 


CMF*FAR 
F*(HT?-HT1 )+FMF*{ HT7-HT6) ) /TMF 
HT4  3 

T43/ETAT 
OHT43I 

K )-(HT3-HT4l ) /CPH 
) -(HT3-HT4) /CPH 
(  (TT4I/TT3(K)  )**(GH/(GH-1.'1)  )  ) 

SG).LT.?.C)  GO  TO  24 
LOG(TT4/TT3(K) )-R*ALOG( PT4/PT3) 


♦  ST3 


GAS  GENERATOR  NOZZLE 


PS5  = 
PS5I 
TT5I 
TS5  = 
TT5  = 
ST5  = 
VA5  = 
A5  =  S 
AMN5 
STAO 
IF(A 
AMN5 
VA5  = 
A5  =  V 
TS5  = 
TT5I 
PS5  = 
ST5  = 
DS5  = 
STAD 


psn 

=  PS 

=  TT 

TT4 

TT4 

CPH 

SQR 

QRT 

=  VA 

0='> 

MN5 

=  1. 

SQR 

A5 

A5* 

=  (F 

PT4 

CPH 

PS5 

0=( 


4*( ( PS5I/PT4)**( (GH-1.0 )/GH)  > 
-  ETAN*<  TT4-TT5I  ) 

*ALOO(TS5/TT3(K > )-R*ALOG< PS5/PT3)  ♦  ST3 

T(2.C*CPH*(TT5-TS5) ) 

(GH*R*TS51 

5/A5 

.0 

.LE.l.C)  GO  TO  15 

0 

T((2.C*CPH*TT5)/(  1.0«-2.0*CPH/(GH*R  )  )  I 

A5/(0F*R) 

TAN*1T4-TT4+TS5)/ETAN 

*( TT5I/TT4)**(GH/(GH-1.( ) ) 

*ALnr,(TS5/TT3lK  )  )-o*ALnG(  PS5/PT^»)     ♦    ST3 

/(R*TS5) 

1.0+FAR  )/( (  l.C+BY)*DS5*VA5)*( PS5-PS0) 


') /2.r *AMN7*AMM7) 
HRFFA 


FAN  MIXING  OLCT 

00    3  1    N=lfNWGR 

MP=WGR (N) /BY 

UMF=WGR (N)*C^F 

TS7=TT7/( l.t  *(GC-1.' 

HS7=( TS7-TSP)*CPC    «■ 

A7=S0RT(GC*R*TS7) 

VA7=AMN7*A7 

VM7=(  VA7  +  MRX-VW7)  /(  1.0  +  MR  ) 

PS7=PT7/(l.C+(GC-l.r) /2.A*AMN7*AMN7)**(GC/{GC-l.O) ) 

DA7=PS7/(R*TS7> 

MIXING    DUCT    TEMPERATURE    ITERATION 

DP     30     IJ  =  l,Nwr>PR 
PSR=PSRR7(  IJ  )*°S7 
0*7=1.9378 
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AN  IV  G  LEVtL  It  MOD  C 


MAIN 


OATF  =  68081 


14/37/^ 


V*8  = 
DO  2 
VA8  = 
VW8  = 
TS8  = 
ZA  = 
ZP  = 

zc  = 

OP    2 
lft    TW8= 
?E  = 

X=(3 

TK  =  ( 

PPV3 

*)*14 

Z    =P 

PFUN 

nze  = 

DT  =  5 
DX  =  - 

ZM  =  { 
*  )  )  /< 

ZN=< 

OPPV 
DZ  =  ( 
9PFU 
TS8  = 
IF(  P 
GO  T 
17  T^3  = 
x=(3 
TK  =  ( 
PPV8 
*)*14 
XP  =  P 


V*7    ♦ 

9     IK  = 

V*8/< 

VR8WA 

530. 

}.243 

.r0  58 

.GOU- 

..^021 

8    IL  = 

TSRBVy 

CPC*< 

MR»VW 

74.11 

273.1 

=       H 

.6959 

PV8/I 

8=ZE/ 

-CPC 

./9.* 

5./Q. 

(  1.     + 

1.  +  7D 

TK*l)X 

8=    PP 

(  P  S  8- 

N8=  (H 

TS8    - 

FIJN8. 

n   16 

TSP8W 
74.  11 
273.1 
=  1C 
.6959 
PV8/  ( 


(PS7 

1 ,  N'wn 

(  l.+M 
{  IK  )* 

7814 
6826 
00C11 
87846 

l  ,nmo 

A(  ID 
TT7-T 
8*V^8 

-  (T 
6  +  ( 
.  **(  A 

♦  144. 

PSP-P 
(HFG( 
-  VR* 
TSRPW 

*  T  S  P  8 
701=  x  ) 
*X  )  ** 

-  X* 
V  8  *  (  X 
PPV8I 
FG(  TW 

PFLN 
LE.C. 


-PS8)*( 1 

VR 

R*VRBWA< 

VA8 


70  2  3  79 

2 

TR 

*TS8 

S8)  +  MP 

l~>.     -     VA 

SK8WA( IL 

TSR8WA( I 

L0G1P(2  1 


./(l.+MR))*(l./(PA7*VA7) +MR/I DW7*VW7) ) 
IK)  )  /(  1.  +MP.)  ) 


*5./9. ) 

)*5./9. ) 

*(<ZA+7B*X+ZC*X**3)/(1.0+7n**)) 


PV8 
TW8 
DHF 
A(  I 
*A( 
*(  Z 
~> 

DTI 
/TK 

*r>p 
8)  * 

8/0 

1E- 


»*n.6 

)  *CJ* 
(  TW8) 
L) 
IL  ) 

b*dx 


*(HF( TW7)-HF( TW3) >*CJ*G  +  MR*V W7*VW7/2 . 
8* V A 8/2. 
)*TS8  -  4Q?. ) 
L)*TS8  -  492. 
8.167)  -  X/TK 

22 

G)  -  Z 

*CJ*G 


+  3.*ZC*OX*X**2)  -  (ZA  +7R*x+  7C*X**M*(  70*DX 


/TK** 
*ZM  + 
PV8  + 
CJ*G* 
PFUM8 
03)  GO.  TO  17 


( ZA+ZR*X+ZC* 

PPV3*DPPV8) / 

0ZE-7F*0HFG(T 


A(  IL  ) 

-  (T 

6  +  ( 

.  **  (  A 
*144. 
PS8-PPV8) *^.6?2 


*TS8 

SR3WA( IL 
TSR8WAI I 
L.nGK<21 


)*TS3  -  4^2. ) 
L)*TS8  -  49?. 
8.167)  -  X/TK 


X**3) /( l.+ZD*X)*7N) 
(  PSB-PPVR  )**2*^.677 
wr  |*(3  )*G)/  (HFG(  TW8)*r.J*O)**?-07 


*5./9.  ) 

)*5./9. ) 

*((7A  +  ZB*X  +  7C*X**3)/<l.'>  +  Zn*X)) 


FAN  NOZZLF 


PS9- 
TS9  = 
DP  ? 
IB  Tl*0  = 
70  =  ( 

*  (S 
X=  (3 
TK  =  ( 
PPV9 

*  )*14 
ZR  =  P 
PFUN 

nzo  = 

PT  =  5 
DX=- 
ZV  =  i 

*) )/( 

OPPV 

nz  =  < 

OPFU 

*G  )  ** 

TS9  = 

IF(P 

GP  T 

19  TW9= 

x=(3 

TK  =  ( 

PPV9 

*)*14 

X9  =  ( 

DP  2 

VA9  = 


PSP 

5  20. 

7  IN= 
TSR9W 
CPC*A 
F(  TW8 
74.  11 
2  73.  1 
=  1C 
.6959 
PV9/( 
9  =  Z0/ 
-CPC/ 
./9.* 

5./q. 
( 1.  + 

l.  +  ZD 
TK*DX 
9=  PP 
(  PS9- 
N9=(  ( 
2-0 
TS9  - 
FIJNQ. 
0  13 
TSR9U 
74.  11 
273.1 
=  1( 
.6959 
PPV9/ 

6  IM  = 
SQRTI 


l.NF 
A(  IN 
LCG( 
)-SF 

-  ( 
6    + 
.**( 
*144 
pcq- 

(  SFG 

TS9 

TSR9 

*TSF 

ZO*x 

*  x )  *■ 

-  X 
VQ*( 
PPV9 
SFGf 
Z 

PFL 
LF.r 


NTR 
)*TS9 
TS8/TS9) 
( TW9)  )*CJ 
TSR9WA( IN 
(TSR9WA( I 
AL0G1C(21 


-    R*ALOG( PS8/DS9)     +    XB*^CG( TW8 )*CJ*G    +     mb* 

*GI 

)*TS9     -     492.) 

M)*TS9     -     492. 

3.167)     -     X/TK 


*5./9.  ) 
)*5. /9. ) 

*((ZA  +  ZR*x  +  zC*x**3)/(l.:-»-7n*x)) 


2/TSR9WA( IN) 

G)     -     ZR 

TW9)  *CJ*G*TSRC»WA(  IN) 


PPV9) *.62: 

( TW9) *CJ*' 

-     MR*DSF( 

WA( IN) 

9WA( IN) 

)*(7B*0X    +    3.*7C*OX*X**7  )     -     (7A    «-ZB*X*7C*X**3  i  *  (  7D*D* 

*2 

*0T  )/TK**2 

x/TK*7M    -f 

)*0PPV9    + 

TW9)*DZ0- 


( ZA+ZB*X+ZC* 

PPV9*DPPV9) / 

Z0*OSFG( TW9)* 


X**^\)  /(  l.  +  ZO*X)*ZN) 
(PS9-PPV9)**2*0.622/TSR9WA(IM) 
TSR9WA( IN) »*CJ*G1 /(SFG(TW9)*CJ* 


N9/0PFIJN9 
.  1E-0  3)    G 


A(  IN)*TS9 

-     (TSR9WA(IN 
6    ♦     (TSR9WA(I 
.**(AL0G1C(21 
*144. 
{  p<$q_ppv9  )*#ft 

1 .NFNVP 
<CPC*(TS8-TS9 


0    TO     19 


)*TS9  -  492. ) 
N)*TS9  -  492. 
8.167)     -     X/TK 


*e5./9.  ) 

)*5./9. ) 

*((ZA  +  ZB*X*ZC*X**3)/(l.';  +  ZD*v)) 

) 


22) /TS"9WA( IN 

)     +    X8*HFG(TWB) *CJ*G    +    M«*HF(TW8)*CJ*G    -    X9* 
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M  IV  G  LFVEL  1,  ^OD  C  MAIN  DATE  =  680*1  14/^7/57 

*  HFG( TH9)*CJ*G  -  MR*HF{  TW9)*CJ*G  +VA8*V  A8/2  .  *  (  I  .  +X8*  (  MR-X8  I  *VR  fi  WA( 
*IK  >**?)  )*2./  (  I.  +  X9  +  ( MR-X9)*VR9WA( IM)**2I > 
VW9  =  VRC>WA(  IM)*VA9 
A9=SQRT(GC*R*TS9) 
AMN9=VA9/A9 
C 
C      SPECIFIC  THRUST  AND  PROPULSIVE  EFFICIENCY  CALCULATIONS 


C 


ST_HRJJS  =  Ii  lL*tFAR  )/(  l.C+BYI)*VA5    -  .1  1.  0/(  1.0*BY!1*QV_ + 
*<BY*BY*X9) /( l.+BY)*VA9    -    (  BY/  (  1  .  +BY  )|  *0V    *■     (  WGR  (  N  )  -BY*X9  )  /  (  1  .  +  RY  )  * 

*  VW9    -     ( WGR(M)  /(  1,+BY)  )*VW0    ♦     STADD)     /G 
THRUST=STHRUS*(CWF    ♦    FMF)     *    G 
TP=THRUST*OV 

ETAPR0=TP/(  TP«-TMF/2.0*(  VA5-OV)**2    +     (  FMF+X9*FMF  )  /2  .*  <  VA9-0V  )  **?    ♦ 

*  ( WMF-FMF*X9)/2.*( VW9- VWO ) **2 ) 
FUMF=r  TMF    -    C^F) 
FUWF=FUMF*G 

SFC=(FUWF /THRUST)* 3600.0 

WRITE(6,20) 
22    FORMAT! //A/, Tft,«  PSR87' » TIB ♦ • VR8WA' , T28 , • TSR 8WA • » T38t • VR9WA* , T4R, • T 
*SP9WA« ,T63,,eR • ,T69, •WGR1 » 

WRITE16, 21 )PSR87( I  J) , VR8WAI IK)tTSRRWAUL) ,VRQWA( IM) f TSR9WAI  IN) ,BY, 
*W6R(N1 

21  FORMAT(2X,7F10.2) 
WPITE(6,22) 

22  FORMAT  (/♦T8»«VA8«fT16,,VWR',T24,«VA9«tT32,«VW9« , TAT t «TSR»fT4R»,T 
*WR',T56,,TSQ«,T64,»TW9«fT71t •STHRUS' ,T81t,AMN'5«,Tfl8,» AMN9' ,  T94,  •  ET 
*APRa*tT10  2»,SFC,,T109,«FAR«,T116,,X9M 

WRITE! 6,2?) VA8tVW8,VA9,VW9,TS8,TWS,TS9,TW9,STHRUS. AMN5 , AMN° . ET APR  0 
*,SFC »FAR,XQ 
2  3  F0RMAT(4X,8Fe.2tF9.2t5F7.3,F7.4) 
GO  TO  26 

24  WRITE(6,25> 

25  FORMAT! T5, 'TURBINE  EXIT  PRESSURE  IS  10WER  THAN  ATMOSPHERIC) 

26  CONTINUE 

27  CONTINUE 

28  CONTINUE 

29  CONTINUE 

30  CONTINUE 

31  CONTINUE 

32  CONTINUE 

33  CONTINUE 

34  CONTINUE 

35  CONTINUE 

36  CONTINUE 
RETURN 
END 
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